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ABSTRACT 

The increasing demand of data, along with the spectrum 
scarcity, are motivating a urgent shift towards exploiting 
new bands. This is the main reason behind identifying 
mm Waves as the key disruptive enabling technology for 5G 
cellular networks. Indeed, utilizing new bands means facing 
new challenges; in this context, they are mainly related to 
the radio propagation, which is shorter in range and more 
sensitive to obstacles. The resulting key aspects that need 
to be taken into account when designing mm Wave cellular 
systems are directionality and link intermittency. The lack 
of network level results motivated this work, which aims at 
providing the first of a kind open source mm Wave frame¬ 
work, based on the network simulator ns-3. The main focus 
of this work is the modeling of customizable channel, phys¬ 
ical (PHY) and medium access control (MAC) layers for 
mm Wave systems. The overall design and architecture of 
the model are discussed in details. Finally, the validity of 
our proposed framework is corroborated through the simu¬ 
lation of a simple scenario. 

Categories and Subject Descriptors 

1.6.5 [Simulation and Modeling]: Model Development— 
Modeling methodolgies; 1.6.7 [Simulation and Modeling]: 
Simulation Support Systems —Environments 

General Terms 

Simulation, Modeling, Architecture, Design, Performance, 
Cellular, Wireless 

Keywords 

mm Wave, 5C, Cellular, Channel, Propagation, PHY, MAC. 

1. INTRODUCTION 

The ever increasing demand of wireless cellular data has mo¬ 
tivated researchers to investigate the potentials of millimeter 
wave communication for the 5th generation of cellular tech¬ 
nology. A substantial body of literature is currently avail¬ 


able discussing physical measurements and formulations for 
millimeter wave channels 0, E As a logical next step, we 
aim at studying how the upper layers of the communication 
network stack work over millimeter wave physical channels. 
In this work we aim to develop the first millimeter wave 
module for the ns-3 network simulator that can be used 
to quantitatively analyze the performance of transport and 
application layer protocols over millimeter wave last-mile 
links. 

The ns-3 network simulator currently implements a wide 
range of network protocols across various layers of the com¬ 
munication network. Due to this it is a valuable tool for 
researchers working on cross-layer design. The ns-3 sim¬ 
ulator already hosts modules for the simulation of WiFi, 
WiMAX and 3CPP-LTE networks. In this paper we pro¬ 
pose the first ns-3 module for the simulation of millimeter 
wave based communication systems. 

The simulation module described in this paper is designed 
to be a fully customizable model where the user can plug 
in various parameters, like carrier frequency, bandwidth, 
frame structure, etc., describing the behavior of the millime¬ 
ter wave channel and devices. In fact, the aim of this work is 
to enable researchers to flexibly use this module for various 
scenarios without the need of altering the source code. 

The rest of the article is organized as follows. In Section 
we discuss the architecture of the mmW module. Sectionl^ 
discusses the modeling of the physical layer. Section fol¬ 
lows with a discussion on the medium access control (MAC) 
layer for our module. The interfacing of the various modules 
is discussed in Section In Section we show the results 
obtained for a simple simulation scenario. Finally, we list 
our future work items and conclude the paper in Section 

2. MMWAVE FRAMEWORK 

Our framework includes a basic implementation of mmWave 
devices, which comprises the propagation and channel model, 
the physical (PHY) layer, and the MAC layer. The module 
completely is developed in C++. The design of this module 
is inspired by the ns-3 Lena module which, in our opinion, 
has a very robust architecture. Fig. gives the UML dia¬ 
gram of the most important classes in our mmWave module. 
We will provide details about the various classes and their 
interoperation in the following sections. 
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Figure 1: Class diagram for the mmWave module. 






























































































































































3. PHY LAYER 

The salient features of the mmWave PHY layer are: (i) a 
fully customizable time division duplex (TDD) frame struc¬ 
ture, (ii) a radio characterization that includes small and 
large scale channel variations, along with supporting multi¬ 
ple input multiple output (MIMO) techniques such as beam¬ 
forming, (iii) a decoding error model at the receiver, (iv) 
an interference model, and (v) a feedback loop for channel 
adaptation. The following parts give a detailed outline of the 
implementation specifics of the PHY layer for the mm Wave 
module. 

3.1 Frame Structure 

The authors in and contend that in order to reduce the 
latency over the air interface, the 5G mm Wave systems will 
be targeted towards TDD operation. The ns-3 module for 
mm Wave implements a customizable TDD frame structure. 
Each frame is subdivided into a number of subframes of 
fixed length specified by the user. Each subframe in turn is 
split into a number of slots of a fixed duration. Each slot 
comprises a specified number of OEDM symbols. A slot can 
be either control or data, assigned for either uplink (UL) 
or downlink (DL). 

Fig. shows an example of the TDD frame structure based 
on the work in [^. Each frame of length 10ms is split in time 
into 10 subframes each of duration 1ms. Each subframe is 
further divided into 8 slots where each slot is of length 125/is 
representing 30 orthogonal frequency division multiplexing 
(OEDM) symbols of length approximately 4.16/is. The first 
two slots are assigned for control in the downlink and in the 
uplink direction respectively. Slots 3 to 5 are allocated for 
downlink and slots 5 to 8 for uplink data transmission. A 
switching gap of 1/is is introduced each time the allocated 
direction changes from uplink to downlink or vice versa. In 
the frequency domain the entire bandwidth of IGHz is di¬ 
vided into 4 resource blocks (RBs). Each RB is subdivided 
into 18 sub-bands each of width 13.89 MHz making a to¬ 
tal of 72 sub-bands for the entire bandwidth. Each of these 
sub-bands is composed of 48 sub-carriers. In this case, the 
total number of resource elements for one slot would be: 
Nre = 30 X (48 X 72) = 103680. 

3.1.1 Parameter configuration 

The frame structure is completely customizable by the user. 

A common object of the mrnwavePhyMacCommon class stores 
the user specified values of all the parameters used by the 
simulator. The parameters used to customize the frame 
structure are specified in Table 

The value of the Transmission Time Interval (TTI), which 
in our model is the duration of one slot, can be derived from 
the parameters in Tableaus, 

TTI = SymbolPerSlot x SymbolLength. 

This is implemented by the function mmwavePhyMacCom- 
mon::GetTTI(). 

Similarly we can compute the bandwidth of one RB as, 

Brb — SubbandsPerRB x SubbandWidth. 



Eigure 2: Example of mmWave frame structure. 


This is implemented by the function mmwavePhyMacCom- 
mon::GetRBWidth(). The total system bandwidth can thus 
be computed as, 

Bsystem = Brr X NumResourccBlock. 

This is implemented by the function mmwavePhyMacGom- 
mon::GetSystemBandwidth(). 

3.1.2 Transmission schemes 

The mmWaveEnbPhy and the mmWaveUePhy models the 
physical layer for the base station and the user device respec¬ 
tively. Broadly the physical layer (i) handles the transmis¬ 
sion and reception of signals, (ii) simulate the start and the 
end of frames, subframes, and slots, (iii) deliver data packets 
and control messages received over the channel to the MAG 
layer, (iv) model the decoding error for the received signal 
and calculate the metrics like the signal to interference and 
noise ratio (SINR). 

The physical layer contains one instance of the mm WaveSpec- 
trumPhy class. The transmission procedure for data frames 
is performed by the StartTxDataFrames. Eor control mes¬ 
sages the StartTxGontrolFrames is invoked. The reception 
procedure is performed by StartRx method. The function¬ 
ality of the StartRx is further subdivided into StartRxData 
and StartRxGontrol. Based on the total band of frequency 
available for transmission, the PHY computes the transmis¬ 
sion power spectral density using the mm WaveSpectrum Val- 
ueHelper component and uses this value for the signal trans¬ 
mission. 

After the reception of the data packets, the PHY layer cal¬ 
culates the SINR of the received signal taking into account 
the MIMO beamforming gains. The physical layer at the 
user device maps the calculated SINR into a Ghannel Qual¬ 
ity Indicator (GQI), which is fed-back to the base station for 
the resource allocation. Gontrol signals are assumed to be 
ideally transmitted. As discussed in Section [33] the PHY 
layer also incorporates the error model where a probabilistic 
approach is used to determine whether a packet should be 
dropped by the receiver. The correctly received packets are 
sent to the MAG layer using the service access point (SAP) 
discussed in Section 0 











































Parameter Name 

Default Value 

Description 

SymbolPerSlot 

30 

Number of OFDM symbols per slot 

SymbolLength 

4.16/is 

Length of one OFDM symbol in /as 

Slots Per Sub frame 

8 

Number of slots in one subframe 

SubframePerFrame 

10 

Number of subframes in one frame 

NumRefereneeSymbols 

6 

The number of reference OFDM symbols per slot 

TDD ControlDataPattern 

“ccdddddd” 

The control (c) and data(d) pattern 

SubearriersPer Subband 

48 

Number of subcarriers in each sub-band 

SubbandsPerRB 

18 

Number of sub-bands in one resource block 

SubbandWidth 

13.89e6 

The width of one sub-band in Hz 

NumResouree Block 

4 

Number of resource blocks in one slot 

C enterFreq 

28e9 

The carrier frequency in Hz 


Table 1: Parameters for configuring the mm Wave frame structure. 


The physical layer handles the start and the end of frames, 
subframes and slots based on the TTI derived from the user 
specified configurations. Based on the resource allocation 
scheme decided by the mmWaveEnhMac^ as described in 
Section [4.2.2[ the PHY decides the nature of the communica¬ 
tion for a particular slot (data/control) and the direction of 
the message to transfer (uplink/downlink). At the beginning 
of each slot the eNodeB PHY sends a Subframeindication to 
the MAC. The subframe indication for the first slot triggers 
the scheduling and resource allocation functions. The data 
packets and the control messages received from the MAC are 
stored in the Packet Burst Queue and the ControlMessage- 
Queue, respectively, and are transmitted to the connected 
device in the allocated slots. 

3.2 Channel Modeling 

As illustrated in Fig. we need to take into account a num- 
ber of procedures to capture characteristics of the mm Wave 
propagation. The key contribution here relates to the com¬ 
putation of the multi-antenna gains, which is particularly 
critical for mm Wave communications. 

The link budget for the mm Wave propagation channel is 
given by, 

+ Gbf -PL- SW, (1) 

where Prx is the total received power in dBm, Ptx is the 
total transmit power, Gbf is the beamforming gain, and 
finally PL and SW represent the pathloss and shadowing, 
respectively. 



Figure 3: mm Wave channel model. 


3.2.1 mmWave Propagation Loss Model 

The mm Wave pathloss model can be modeled with 3 states, 
as reported in [^, line of sight (LOS), non-line of sight (NLOS) 
and outage. For each link, the channel is determined through 
the following procedure: 

• based on the distance between the UE and the eNB, de¬ 
termine the probability of the link being in each of the 
three states using the model in {PloS, PnloS 5 Pout)'l 

• uniformly pick a reference value (Pref) between 0 and 
1 and compare with the probability associated with 
each channel state; 

• if Pref < Pros, pick LOS channel; 

else if Plos < Pref < Pros PPnlos, pick NLOS 
channel; 

otherwise pick outage. 

For each link, on determining the channel state, the pathloss 
and shadowing is obtained by, 

PL{d)[dB]=a +I310logio{d) + ^, (2) 

where ^ represent shadowing, parameter d represents the dis¬ 
tance from receiver to transmitter, the values of parameter 
a, /3, and a for each channel scenario are given in [^. 

3.2.2 MIMO 

Due to the high pathloss, multiple antenna with beam form¬ 
ing is essential to ensure acceptable range of communication 
in mm Wave systems. We briefly discuss the associated con¬ 
cepts of channel matrices and beam forming in this respect. 

Channel matrix: We model the mm Wave channel as a 
combination of clusters, each composed of several subpaths Q 
The channel matrix is described by the following: 

HitJ) = ( 3 ) 

k=l1=1 

where, K is the number of clusters, Lk the number of sub¬ 
paths in cluster k, gki (t, /) the small-scale fading over fre¬ 
quency and time; Urx(')is the spatial signature of the re¬ 
ceiver, Utx(-) the spatial signature of the transmitter. 

^See for a long term statistical characterization of the 
mmW^e channel. 































































The small-scale fading is generated based on the number 
of clusters, number of subpaths per cluster, Doppler shift, 
power spread, delay spread and angle of arrival (AoA) as 
given in 1 by, 

9H{t, f) = ( 4 ) 

where, Pki is the power spread; fd is the maximum Doppler 
shift; (jJki is the AoA of the subpath relative to the direction 
of motion; Tki gives the delay spread, and / is the carrier 
frequency. 

Beamforming: In order to support phased-array antennas, 
a new Antenna Array Model class is developed, which con¬ 
tains a complex beamforming vector. For both transmitter 
and receiver, based on the channel matrices, the beamform¬ 
ing vectors are computed using the power algorithm. 

The beamforming gain from transmitter i to receiver j is 
given as, 

G{t, f)ij = |w;„y H(t, 1^ (5) 

where, H(i, f)ij is the channel matrix of link, is 

the beamforming vector of transmitter i, when transmitting 
to receiver j and Wrxij the beamforming vector of receiver 
j, when receiving from transmitter i. 

3.2.3 Channel Configuration 

To reduce the computational complexity, the channel matri¬ 
ces and beamforming vectors are pre-generated in MATLAB®. 

Load Files: At the beginning of each simulation we load 
100 instances of the spatial signature matrices, along with 
the beamforming vectors. Then, a channel matrix instance 
per UE-eNodeB pair is randomly picked to characterize the 
radio link. As we will discuss later, we simulate the long 
term fading by randomly picking an instance of the channel 
form the pre-generated files. 

Link Initialization and Channel Matrix Updates: In 

the mmWaveBeamForming class we define the a member 
m_ehannelMatrixMap (Fig. to map the channel ma¬ 
trix instance to each radio link. During the simulation, the 
small-scale fading is calculated at every slot, based on Eq|^ 
The speed of the user is obtained directly from the mobility 
model. The remaining parameters are only subject to the 
environmental conditions , rural, urban etc., and therefore 
assumed constant over the entire simulation time. The refer¬ 
ence values for different environments were recorded during 
the our mm Wave channel propagation measurement cam¬ 
paign reported in [^. On the other hand, for the large-scale 
fading, the spatial signature matrices are periodically up¬ 
dated with customizable interval, say 100 ms, to simulate a 
sudden change of the perceived channel. 

Beamforming Vector: When configuring each radio link, 
the beamforming vector is stored in the antenna array model 
of both base station and user side; the former will store the 
beamforming vectors of all the UEs, the latter will only store 
the beamforming vectors of the base stations within range. 

The beamforming parameters are defined using following 
structures: 
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Eigure 4: MIMO channel gain over frequency and 
time. 


struct BeamEormingParams 

{ 

complexVector_t m_txW; 
complexVector_t m_rxW; 

ChannelMatrix m_channelMatrix; 

} 

struct ChannelMatrix 

{ 

complex2DVector_t m_txSpatialMatrix ; 
complex2DVector_t m_rxSpatialMatrix ; 

} 


We can observe in Fig. [^the channel gain trend obtained 
in a scenario where the number of antennas at the base sta¬ 
tion and mobile device is 64 and 16, respectively, and the 
user is moving at a speed of 36 km/h. The beamforming 
gain with small-scale fading is calculated using on Eq. 
The small-scale fading includes two components, frequency 
selective fading and time selective fading, caused by multi- 
path effect and Doppler effect, respectively. The spatial sig¬ 
natures and beamforming vectors are periodically updated 
with larger interval to capture the effects of long term fading 
which cause the sudden drop of beamforming gain around 
100 ms. 


3.3 Error Model 

Similar to ns-3 Lena module, our mm Wave module includes 
a error model for data packets according to the standard 
link-to-system mapping (LSM) techniques. By utilizing the 
LSM and the Mutual Information Based Effective SINK 
(MIESM) [^, the receiver computes the error probability 
for each transport block (TB) and determines whether the 
packet can be decoded or not. The TB can be composed of 
multiple codeblocks (CB) and its size depends on the chan¬ 
nel capacity. The block error probability (BLER) of each 
CB depends on its size and associated MCS: 


CBLER,i{'yi) — 


1 


\ — 0j-f f 7^ bcsizE^MCS \ 

V ^s/^C-C!QJ2;E .fMCS J 


(6) 


where 7 ^ is the mean mutual information per coded bit 
(MMIB) of the codeblock i, hcsizEMCS and ccsizeMCS 
corresponds to the mean and standard deviation of the Cans- 







Figure 5: Interference model. 


sian cumulative distribution, respectively. Now we can com¬ 
pute the TB block error rate: 

c 

TbLER = 1 — ]^(1 “ CBLER.ii'Ji))- ( 7 ) 

i=l 

In case of failure, the PHY layer does not forward the in¬ 
coming packet to the upper layers and, at the same time, 
triggers a retransmission process]^ 

3.4 Interference 

Albeit being presumably less threatening in the mm Wave 
regime, because of the directionality of the multiantenna 
propagation, interference computation is still pretty relevant 
in terms of system level simulations. In fact, there might be 
some special spatial cases where interference is non negli¬ 
gible. Therefore, we propose an interference computation 
scheme that takes into account the beamforming directions 
associated with each link. 

We will use Fig. m as a reference. As an example, we 
compute the SINK between BSi and UEi. To do so, we 
hrst need to obtain the channel gains associated with both 
the desired and interfering signals. Following Eq. we get 

Gii = 

G 2 I = f)21-Wtx22f ■ 


We can now compute the SINK: 


SINRii = 


Ptx,11 

PLii 


Gii 


Ptx,22 

PL 2 I 


G 2 I 


+ BW X No ’ 


( 9 ) 


where Ptx,ii is the transmit power of BSi, PLu is the 
pathloss between between BSi and UEi^ and BW x Nq is 
the thermal noise. 


3.5 CQI Feedbacks 

In order to ensure reliable communication over a variable 
channel, feedback mechanisms are key in mostly all cellular 
communication systems. Similar to LTE we utilize the CQI 
feedback scheme for our module. The downlink CQI feed¬ 
back message is generated by the mmWaveUePhy. In our 
simulator the UE computes the CQI based on the SINK of 
the signal received in a particular data slot. The computa¬ 
tion of CQI is the same as that for LTE as given in [ 9 ] and 

m _ 

^This can be a TCP retransmission, or an hybrid automatic 
repeat request (HARQ), which is part of our future work. 


4. MAC LAYER 

The MAC layer is developed using the class mmWaveMac 
which is the base class for the mmWaveEnhMac for the eN- 
odeB and the mmWaveUeMac for the user. The chief func¬ 
tion of this layer is to deliver data packets coming from the 
upper layers (the net device in this case) to the physical 
layer and vice-versa. In fact this layer is designed for the 
synchronous delivery of upper layer data packets to the PHY 
layer which is key for proper data transfer in TDD mode. 

The eNodeB MAC layer is connected to the scheduler mod¬ 
ule using the MAC-SCHED service access point (Sec. |5.2| . 
The relationship between the PHY, MAC and the scheduler 
module for the cNodeB is depicted in Eig. Thus the 
MAC layer communicates the scheduling and the resource 
allocation decision to the PHY layer. The scheduler hosts 
the adaptive modulation and coding (AMC) module. The 
following sub-sections discuss these features in depth. 

4.1 Adaptive Modulation and Coding 

The working of the AMC is similar to that for LTE. The user 
measures the CQI for each downlink data slots it is allocated. 
The CQI information is then forwarded to the eNodeB us¬ 
ing the mmWaveCqiReport control message. The eNodeB 
scheduler uses this information to compute the most suit¬ 
able modulation and coding scheme for the communication 
link. 

The AMC is implemented by the eNodeB MAC schedulers. 
During resource allocation, for the current frame work, the 
wide band CQI is used to generate the modulation and cod¬ 
ing scheme (MCS) to be used and the transport block (TB) 
size that can be transmitted over the physical layer. The 
AMC module provides this frame work for the unique map¬ 
ping of the CQI, the MCS, spectral efficiency and the TB 
size. The TB size is calculated based on the values of the 
total number of subcarriers per resource block derived from 
the user customized conhguration, the number of symbols 
per slot and the number of reference symbols per slot. A 
cyclic redundancy code (CRC) length of 24 bits is used. 

4.2 Scheduler 

Eollowing the design strategy for the ns-3 LTE module pR] , 
the virtual class mm WaveMacScheduler dehnes the interface 
for the implementation of MAC scheduling techniques. The 
scheduler performs the scheduling and resource allocation 
for a subframe with both downlink and uplink slots. 

4.2.1 TDD scheme 

The TDD scheme enforced by the scheduler module is based 
on the user specihed parameter “TDDControlDataPattern” 
given in Table The slots specihed for control are assigned 
alternately for downlink and uplink control channels. The 
data slots are equally divided between downlink and uplink 
slots with the hrst n/2 data slots allocated to downlink data 
and rest for uplink, where n is the total number of data slots. 
This scheme minimizes the switching time between uplink 
and downlink data transmissions. 

The scheme described in Eig. [^is an example of the imple¬ 
mentation of the above algorithm with the default control 
data pattern. This module will be further enhanced in future 
to incorporate more advanced features like dynamic TDD. 












Figure 6: PHY, MAC and scheduler modules with 
the associated SAPs. 

4.2.2 Resource Allocation 

Using the mmWavePhyMacCommon object, the division of 
resources in the frequency domain can be customized as 
given in Table 

The MAC scheduler currently implements a simple round 
robin algorithm to allocate uplink and downlink data slots 
to the connected users. All the frequency elements in a 
particular slot are assigned to the same user. The control 
slots are not allocated to any particular user. Any or all 
user can receive from and transmit to the base station in 
the control slots. 

For the case of the round robin scheduler the use of the 
CQI is limited to the determination of the transport block 
size. In future the CQI information will be used to actively 
control the scheduling decisions. 

5. SERVICE ACCESS POINTS 

The interface between the PHY and the MAC and the MAC 
and the scheduler are dehned as service access points (SAPs) 
as given in [^. The relationship between the PHY, MAC 
and the scheduler through their associated SAPs are shown 
in Fig. The relationship between modules connected 
through SAPs can be viewed as that of service providers 
and users. The SAP provider caters to the requirement of 
the SAP user based on certain requests received from the 
user. 


5.1 PHY-MAC 

The communication between the MAC and the PHY layer 

using the MAC-PHY SAP is through the following processes: 

i Subframe Indication: The subframe indication is sent 
by the PHY layer to the MAC at the beginning of each 
slot (unlike LTE where it is sent every subframe). The 
subframe indication for slot 1 for a particular subframe 
triggers the scheduling procedure for the eNodeB MAC. 
The subsequent indications are required for proper deliv¬ 
ery of upper layer data. 

ii Data transmission: The eNodeB MAC maintains data 
queues for each of the connected UE and just one such 
queue is sufficient for the user device. Based on the schedul¬ 
ing scheme and the allocated resources, the MAC layer 
will send the scheduled number of packets (given by the 
transport block size) to the PHY layer for transmission 
over the radio link. 


iii Scheduling and allocation notifications: The schedul¬ 
ing and resource allocation decision received by the eN¬ 
odeB MAC from the scheduler is relayed to the PHY 
layer using the mmWaveResourceAllocation message. The 
PHY of the base station in turn transmits this message to 
all the connected users notifying all the attached devices 
of the scheduling decision. 

iv CQI notification: Based on the SINK of the received 
data slots, the UE PHY calculate the CQI and transmits 
it to the base station in the next uplink control slot. The 
eNodeB PHY on receiving the mmWaveCqiReport control 
message, relays it to the MAC. 

5.2 MAC-SCHED 

The eNodeB MAC uses the service provided by the scheduler 

by the following processes: 

i Trigger Request and Configuration Indication: On 

receiving the Subframe Indication for slot 1 of a particu¬ 
lar subframe, the MAC sends a Scheduling Trigger Request 
to the scheduler for the (subframeNum + delayY^ sub- 
frame, where the delay is specified by the user using the 
parameter LlL2ControlLatency. The scheduler returns 
the scheduling and allocation decisions in the Scheduling 
Configuration Indication in response to the trigger. 

ii CQI notification: The eNodeB MAC on receiving the 
CQI information from the PHY, sends it to the scheduler. 
The scheduler needs this information for future scheduling 
decisions. 


6. SIMULATION RESULTS 

We validate our framework through a simple simulation sce¬ 
nario, where we consider one eNodeB and one UE. The user 
is at an initial distance of 40m from the base station and 
moves towards the cell-edge with a constant velocity of 20 
m/s. The simulation configuration settings correspond to 
the default values shown in Table On top of that, the 
eNB has a transmit power of 30 dBm and a receiver noise 
figure of 5dB. Simulations are executed twice to capture the 
difference between the line of sight (LoS) and the non line 
of sight (NLoS) scenario. 

The simulation results are reported in Fig. . On the one 
hand. Figs. |7a| and show the variation of the downlink 
SINK with both time and frequency for the LoS and NLoS 
case, respectively. On the other hand, in Figs. |7d| and 
[7c] we can observe the decrease in downlink data rate and 
average downlink SINK with the increasing distance between 
the UE and the base station, for both LoS and NLoS links. 

7. CONCLUSIONS & FUTURE WORK 

In this paper, a novel module for the simulation of mm Wave 
cellular systems has been presented. The module, which is 
publicly available at [^, is highly customizable to facilitate 
researchers to use it flexibly and analyze different scenar¬ 
ios using varying configurations. A basic implementation of 
mm Wave devices, MAC layer, PHY layer and channel mod¬ 
els are developed. As part of our future work, we aim at in¬ 
troducing several enhancements, which include the integra¬ 
tion of (i) higher layer modules for end-to-end performance 
analysis, (ii) a HARQ module, (iii) uplink power control. 



















(a) SINK for user with LoS link. 


(b) SINR for user with NLoS link. 



(c) Average SINR estimated by the users. 


(d) Data rate for users with LoS and NLoS links. 


Figure 7: Simulation results 


(iv) ns-3 based channel matrix generation and beamforming 
computation, (v) more sophisticated MAC algorithms, (vi) 
relay devices and the evaluation of (i) multiple access strate¬ 
gies and (ii) Transmission Control Protocol (TCP) perfor¬ 
mance over mm Wave cellular networks. 
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